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We study the evolution of a flat Friedmann-Robertson Walker Universe, filled with a bulk vis¬ 
cous cosmological fluid, in the presence of variable gravitational and cosmological constants. The 
dimensional analysis of the model suggest a proportionality between the bulk viscous pressure of 
the dissipative fluid and the energy density. With the use of this assumption and with the choice of 
the standard equations of state for the bulk viscosity coefficient, temperature and relaxation time, 
the general solution of the field equations can be obtained, with all physical parameters having a 
power-law time dependence. The symmetry analysis of this model, performed by using Lie group 
techniques, confirms the unicity of the solution for this functional form of the bulk viscous pressure. 
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I. INTRODUCTION 

In many cosmological and astrophysical situations an idealised fluid model of the matter is inappropriate. Such 
possible situations are the relativistic transport of photons, mixtures of cosmic elementary particles, the evolution of 
cosmic strings due to their interaction with each other and with the surrounding matter, the classical description of 
the (quantum) particle production phase, interaction between matter and radiation, quark and gluon plasma viscosity, 
different components of dark matter etc. JjJ. From a physical point of view the inclusion of dissipative terms in the 
energy momentum tensor of the cosmological fluid seems to be the best-motivated generalization of the matter term 
of the gravitational field equations. 

The theories of dissipation in Eckart-Landau formulation , who made the first attempt at creating a relativistic 

theory of viscosity, are now known to be pathological in several respects. Regardless of the choice of equation of state, 
all equilibrium states in these theories are unstable. In addition, as shown by Israel [|J, signals may be propagated 
through the fluid at velocities exceeding the speed of light. These problems arise due to the first order nature 
of the theory since it considers only first-order deviations from the equilibrium, leading to parabolic differential 
equations, hence to infinite speeds of propagation for heat flow and viscosity, in contradiction with the principle of 
causality. While such paradoxes appear particularly glaring in relativistic theory, infinite propagation speeds already 
constitutes a difficulty at the classical level, since one does not expect thermal disturbances to be carried faster than 
some (suitably defined) mean molecular speed. Conventional theory is thus applicable only to phenomena which are 
“quasi-stationary” i.e. slowly varying on space and time scales characterized by mean free path and mean collision 
time ||]. This is inadequate for many phenomena in high-energy astrophysics and relativistic cosmology involving 
steep gradients or rapid variations. These deficiencies can be traced to the fact that the conventional theories (both 
classical and relativistic) make overly restrictive hypothesis concerning the relation between the fluxes and densities 
of entropy, energy and particle number. 

A relativistic second-order theory was found by Israel 0 and developed by Israel and Stewart || into what is called 
‘transient’ or ‘extended’ irreversible thermodynamics. In this model deviations from equilibrium (bulk stress, heat 
flow and shear stress) are treated as independent dynamical variables leading to a total of 14 dynamical fluid variables 
to be determined. The solutions of the full causal theory are well behaved for all times. Hence the advantages of the 
causal theories are the followings ||: 1) for stable fluid configurations the dissipative signals propagate causally 2) 
unlike Eckart-type’s theories, there is no generic short-wavelength secular instability in causal theories and 3) even 
for rotating fluids, the perturbations have a well-posed initial value problem. Therefore, the best currently available 
theory for analyzing dissipative processes in the Universe is the full Israel-Stewart causal thermodynamics. 
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Due to the complicated nonlinear character of the evolution equations, very few exact cosmological solutions of the 
gravitational field equations are known in the framework of the full causal theory. For a homogeneous Universe filled 
with a full causal viscous fluid source obeying the relation £ ~ p *, exact general solutions of the field equations have 
been obtained in [Q, ||, [jllj-jl5). In this case the evolution of the bulk viscous cosmological model can be reduced 
to a Painleve-Ince type differential equation. It has also been proposed that causal bulk viscous thermodynamics can 
model on a phenomenological level matter creation in the early Universe [||, fl(| . 

Recent observations of type la supernovae with redshift up to about z < 1 provided evidence that we may live in 
a low mass-density Universe, with the contribution of the non-relativistic matter (baryonic plus dark) to the total 
energy density of the Universe of the order of ~ 0.3 The value of Q m is significantly less than unity |2(J 

and consequently either the Universe is open or there is some additional energy density p sufficient to reach the value 
fltotal = 1, predicted by inflationary theory. Observations also show that the deceleration parameter of the Universe 
q is in the range — 1 < q < 0, and the present-day Universe undergoes an accelerated expansionary evolution. 

Several physical models have been proposed to give a consistent physical interpretation to these observational facts. 
One candidate, and maybe the most convincing one for the missing energy is vacuum energy density or cosmological 
constant A juj. 

Since the pioneering work of Dirac |2^| , who proposed, motivated by the occurrence of large numbers in Universe, a 
theory with a time variable gravitational coupling constant G, cosmological models with variable G and nonvanishing 
and variable cosmological term have been intensively investigated in the physical literature |^3)-|42 . In the isotropic 
cosmological model of Chen and Wu (30) it is supposed, in the spirit of quantum cosmology, that the effective 
cosmological constant A varies as a~ 2 (with a the scale function). In the cosmological model of Lima and Maia |4d[] 
the cosmological constant A = A ( H ) = 3 /3H 2 + 3 (1 — f3) H 3 /Hi is a complicated function of the Hubble parameter 
H , a constant (3 and an arbitrary time scale HJ , leading to a cosmic history beginning from an instability of the de 
Sitter space-time. The cosmological implications of a time dependence of the cosmological of the form A ~ t ~ 2 have 
been considered by Berman ^ 8 | . Waga (42| investigated flat cosmological models with the cosmological term varying 
as A = a/a 2 + /3H 2 + 7 , with a, [3 and 7 constants. In this model exact expressions for observable quantities can 
be obtained. Nucleosynthesis in decaying-vacuum cosmological models based on the Chen-Wu ansatz |j(J has been 
investigated by Abdel-Rahman p4| . The consistency with the observed helium abundance and baryon asymmetry 
allows a maximum vacuum energy close to the radiation energy today. Anisotropic Bianchi type I cosmological 
models with variable G and A have been analyzed by Beesham p7| and it was shown that in this case there are no 
classical inflationary solutions of pure exponential form. Cosmological models with the gravitational and cosmological 
constants generalized as coupling scalars and with G ~ a n have been discussed by Sistero [Q. Generalized field 
equations with time dependent G and A have been proposed in |k| and J3(| in an attempt to reconcile the large 
number hypothesis with Einstein’s theory of gravitation. Limits on the variability of G using binary-pulsar data have 
been obtained by Damour, Gibbons and Taylor |n[]. A detailed analysis of Friedmann-Robertson-Walker (FRW) 
Universes in a wide range of scalar-tensor theories of gravity has been performed by Barrow and Parsons j26j. 

It is the purpose of the present paper to consider the evolution of a causal bulk viscous fluid filled flat FRW type 
Universe, by assuming the standard equations of state for the bulk viscosity coefficient, temperature and relaxation 
time and in the presence of variable gravitational and cosmological constants. In order to obtain some very general 
properties of this cosmological model with variable constants we shall adopt a method based on the studies of the 
symmetries of the field equations. As a first step we shall study the field equations from dimensional point of view. 
The dimensional method provides general relations between physical quantities and allow to make some definite 
assumptions on the behavior of thermodynamical quantities. In particular we find that, under the assumption of 
the conservation of the total energy of the Universe, the bulk viscous pressure of the cosmological fluid must be 
proportional to the energy density of the matter component. With the use of this assumption the gravitational field 
equations can be integrated exactly, leading to a general solution in which all thermodynamical quantities have a 
power-law time dependence. 

As we have been able to find a solution through dimensional analysis, it is possible that there are other symmetries 
of the model, since dimensional analysis is a reminiscent of scaling symmetries, which obviously are not the most 
general form of symmetries. Hence we shall study the model through the method of Lie group symmetries, showing 
that under the assumed hypotheses of the proportionality of bulk viscous pressure to the energy density, there are no 
other solutions of the field equations. 

The present paper is organized as follows. The field equations are written down in Section II. The dimensional 
analysis of the model is performed in Section III. The general solution of the gravitational field equations for the bulk 
viscous pressure proportional to the energy density is obtained in Section IV. The Lie group symmetry study of the 
model is considered in Section V. In Section VI we discuss and conclude our results. 
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II. FIELD EQUATIONS, THERMODYNAMICS AND CONSEQUENCES 

In the presence of a time variable gravitational and cosmological constants G and A the Einstein gravitational field 
equations are: 


1 „ 8nG(t) 

Rik ndikR — 4 Tik A (t)gik- 


Cl) 


Applying the covariance divergence to the second member of equation ( 0 ) one obtain the generalized “conservation 
equation” 


TL = - 


Gj 

G 


t? - 


c 4 A,i 

8t tG 


( 2 ) 


We assume that the geometry of the space-time is of FRW type, with a line element 

ds 2 = c 2 dt 2 - a 2 it) ( dx 2 + dy 2 + dz 2 ) . (3) 

The Hubble parameter associated with the line element (j|) is defined as H = £. 

The energy-momentum tensor of the bulk viscous cosmological fluid filling the very early Universe is given by 

T? = (pc 2 +P + n) Ul u k - (p + H) 6 k , (4) 

where p is the mass density, p the thermodynamic pressure, n the bulk viscous pressure and iq the four velocity 
satisfying the condition u t u l = 1. The particle and entropy fluxes are defined according to N 1 = nu l and S l = 
eN l — u l , with n is the number density, e the specific entropy, T > 0 the temperature, £ the bulk viscosity 

coefficient and r > 0 the relaxation coefficient for transient bulk viscous effect (i.e. the relaxation time). 

The evolution of the cosmological fluid is subject to the dynamical laws of particle number conservation JVU = 0 

and Gibb’s equation Tde = d +pd (^). In the following we shall also suppose that the energy-momentum tensor 

of the cosmological fluid is conserved, that is T k k = 0. 

For the evolution of the bulk viscous pressure we adopt the causal evolution equation , obtained in the simplest 
way (linear in n) to satisfy the H-theorem (i.e. for the entropy production to be non-negative, S'^ = > 0 !)• 

Therefore the gravitational field equations describing the cosmological evolution of a causal bulk viscous fluid in 
the presence of variable gravitational and cosmological constants are 


3 U 2 = 8irG(t)p + A(t)c 2 , (5) 

2 H + 3 H 2 = - 87rC {p + n) + A (t)c 2 , (6) 

P + 3 (^p+ H = 0, (7) 

8npG(t) + A(t)c 2 = 0, (8) 

rn + n = -3^-|rnf3F + ^-|-^j . (9) 


In eq. ( 0 ), e = 0 gives the truncated theory (the truncated theory implies a drastic condition on the temperature), 
while e = 1 gives the full theory. The non-causal theory has r = 0. Dot denotes differentiation with respect to time. 

In order to close the system of equations (|[]|) we have to give the equation of state for p and specify T, £ and r. 
As usual, we assume the following phenomenological laws Q : 


p = ( 7 -l)pc 2 , £ = ap s , T = f3p r , r = ip 1 = ap s \ 


( 10 ) 


where a > 0, [3 > 0 are dimensional constants and 1<7<2, s>0 and r > 0 are numerical constants. Eqs. (|To[) are 
standard in cosmological models whereas the equation for r is a simple procedure to ensure that the speed of viscous 
pulses does not exceed the speed of fight. These are without sufficient thermodynamical motivation, but in absence 
of better alternatives we shall follow the practice adopting them in the hope that they will at least provide indication 
of the range of possibilities. The temperature law is the simplest law guaranteeing positive heat capacity. 
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With the use of equations of state (@ the evolution equation of the bulk viscous pressure becomes 

P 


n 


V~ s n = - 3 p h - in 

a 2 


3H — (1 + r) 


P 


( 11 ) 


In the context of irreversible thermodynamics p 1 p, T and the particle number density n are equilibrium magnitudes 
which are related by equations of state of the form p = p(T, n) and p = p(T , n). From the requirement that the entropy 

is a state function, we obtain the equation = P ^ C n +P ~ ^ [-§p ( p/c 2 )] n . For the equations of state ([!(]) this 

relation imposes the constraint r = , so that 0 < r < 1/2 for 1 < 7 < 2, a range of values which is usually 

considered in the physical literature jTj]. 

The growth of the total comoving entropy £ over a proper time interval {to,t) is given by (0: 

q /■* TTHn 3 

£(i)-E(io) = -^jf (12) 


where ks is the Boltzmann’s constant. 

The Israel-Stewart-Hiscock theory is derived under the assumption that the thermodynamical state of the fluid is 
close to equilibrium, that is the non-equilibrium bulk viscous pressure should be small when compared to the local 
equilibrium pressure |II| << p = (7 — 1) pc 2 Q|. If this condition is violated then one is effectively assuming that 
the linear theory holds also in the nonlinear regime far from equilibrium. For a fluid description of the matter, the 
condition ought to be satisfied. 

An important observational quantity is the deceleration parameter 9 = 7 ^ (jy) ~ 1- The sign of the deceleration 
parameter indicates whether the model inflates or not. The positive sign of q corresponds to “standard” decelerating 
models whereas the negative sign indicates inflation. 

By using the field equations we can express the gravitational constant in the form 

= —■ (13) 

47r p 

With the use of this expression for G we obtain the cosmological constant in the form 


A(f)c 2 = 3 H 2 



(14) 


From the field equations we obtain for the derivative of the Hubble function the alternative expression 

47rG(f) 2 


H = — 


(, oc 2 + p + n). 


(15) 


III. DIMENSIONAL ANALYSIS OF FIELD EQUATIONS 

In this section we show how writing the field equations in a dimensionless way a particular solution of the field 
equations, corresponding to a specific choice of the bulk viscous pressure, can be obtained. 

The 7 r — monomia is the main object in dimensional analysis. It may be defined as the product of quantities which 
are invariant under change of fundamental units. 7 r — monomia are dimensionless quantities, their dimensions are 
equal to unity. The dimensional analysis has structure of Lie group |45| ]. The 7 r — monomia are invariant under the 
action of the similarity group. On the other hand we must mention that the similarity group is only a special class 
of the mother group of all symmetries that can be obtained using the Lie method. For this reason when one uses 
dimensional analysis only one of the possible solutions to the problem is obtained. 

The equations (|§||) and the equation of state ( 0 ) can be expressed in a dimensionless way by the following 
7 T — monomia: 


7 Tl = 


Gpt 2 


7T 7 


( 

JL 

ap s 


7T2 = 


GW 2 ^ _ r .2 ^ _ |n| 

9 5 ^3 P" 1 7T4 


P 


_ r _ 

7rB_ init’ ^ t ~ rH ' 


(16) 


£ T pc 2 , 2 

^8 = —, 719 = -5-r) Mo = -) Mi = A c t . 

rp (jp r p 


(17) 
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The following relation can be obtained from the it — monomia. From 7 Ti , tt 2 , 7 t 3 , 774 and 7Tio we can see that 


p oc p oc |n|, 


( 18 ) 


note that [|II|] = [p] = [ pc 2 ] , that is, they have the same dimensional equation and therefore from this point of view 
they must have the same behavior in order of magnitude. Now, from 777 = and 7 r s = — we obtain 


implying 


From 7 t 3 and 7 r 8 we find 


or 


^8 


ap 



Gpt 2 ap s 1 


(19) 


( 20 ) 


( 21 ) 


If s = 1/2 we obtain the relationship a 2 = c 2 /G. If s / 1/2 the “constant” G must vary. 
Therefore the solutions that Dimensional Analysis suggests us are the following ones: 

p cx (a _ 1 f) “/pot (a _1 t) 3-1 = (7 — 1) pc 2 , |II| oc (cU 1 !) S_1 = xpc 2 
A oc c~ 2 t~ 2 , Goca^c 2 !”, s/1, 


( 22 ) 


(23) 

(24) 


with x a numerical constant. From the equation of state of bulk viscosity coefficient, temperature and relaxation time 
we find 

£ oc a (a _ 1 t ) S_1 , T oc (3 S_1 , r oc t, s/ 1. (25) 


IV. GENERAL SOLUTION OF THE GRAVITATIONAL FIELD EQUATIONS WITH BULK VISCOUS 
PRESSURE PROPORTIONAL TO THE ENERGY DENSITY 

As we have seen in the previous Section, dimensional analysis suggests us that |II| oc pc 2 . Hence generally the 
bulk viscous pressure can be represented in the form n = — xpc 2 , with x > 0 and x ^ 1- With this hypothesis the 
conservation equation becomes 


P + 3 (7 - X) P H = 0, 

leading to the following relationship between the energy density p and the scale factor a: 

Po 


P = 


,3(7—x) ’ 


(26) 


(27) 


where po > 0 is an arbitrary constant of integration. With the use of the assumed functional form of n the evolution 
equation of the bulk viscous pressure becomes 


p = -K 0 p 2 s , 


where we denoted Kn = — 

u a. 2 


2-x 

2x(T“X) 


The solution of this equation is 


Pit) 


Po 

1 5 


(28) 


(29) 
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with po = [(1 — s ) Kq] 1 ^ 11 1 ^. From Eqs. (27) and we obtain the time dependence of the scale factor 

a(t) = ao t 3 ( 1 -»KT-x). (30) 


with a 0 = (&) • 

The Hubble parameter is 


H(t) 


Ho 
t ’ 


(31) 


where Hq = > 0- The behavior of the gravitational and cosmological constant can be obtained from Eqs. 

© and Jl4|), respectively: 


Git) 


G 0 

1-2 s 


W) = % 


(32) 


where Gq = ^ H - 1 * s 1 and Ao = 3T7 2 (2s — 1 ) /c 2 . The behavior of the bulk viscosity coefficient, temperature and 
relaxation time is given by 


£ = 



T = 


n_ 

t ^ 


T = Tot, 


(33) 


where we denoted £o = app, Tq = Ppq and r 0 = £ 0 /po■ The expression of r is in agreement with the one expected 
from a theoretical point of view, as argued in |43| , since for a viscous expansion to be non-thermalizing we should have 
t < t, for otherwise the basic interaction rate for viscous effects could be sufficiently rapid to restore the equilibrium 
as the fluid expands. 

The deceleration parameter q behaves as: 


q = 1/H 0 — 1 = const. (34) 

For values of the parameters s, 7 and x so that 0 < Hq < 1 the expansion of the Universe is non-inflationary, while 
for Hq > 1 we obtain an inflationary behavior. 

The comoving entropy in this model is 


E(f)-E(t 0 ) = E 0 f^(^+ 1 - 1 ), 


where we denoted Eo = ^ 7 

kBT 0 { — +rUl) 

The ratio of the bulk viscous and thermodynamic pressure is given by 


n 


X 

p 


7-1 


(35) 


(36) 


Since x is assumed to be a small number, the present model is thermodynamically consistent for the matter equations 
of state of cosmological interest, describing high-density cosmological fluids, when bulk viscous dissipative effects are 
important. 

We would like to point out that we have obtained the same results as in the previous section, at least in order of 
magnitude. By direct integration of the field equations we have been able of finding some of the numerical constants 
describing the evolution of the Universe. This is a crucial issue at least for the cosmological constant. With this 
complete solution it is observed that we have the special case s = 1/2 for which we obtain 

G = const., A = 0, poet -2 , a oc t 3 ^-x ),..., E(£) — E {to) oc t 2 G-x +r ~ 1 '>. (37) 


If for example we take 7 = 4/3 (radiation predominance) then 

E(f) — E (f 0 ) oc D. 


(38) 


We therefore see that the x~~ parameter, i.e. the causal bulk effect, perturb weakly the FRW perfect fluid solution. If 
X = 0 we recover the perfect fluid case, and then the comoving entropy is constant. 
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If s ^ 1/2 then we obtain the following behavior for G and A: ifl/2<s<l then G is a growing function on time 
while A > 0, and A oc t~ 2 , if s < 1/2 and A < 0 or s > 1 and A < 0, then G is a decreasing function on time while 
A oc t~ 2 . 

These results agree with the ones obtained in a recent work pd)] where we studied through different routes, renor¬ 
malization group (RG), dimensional analysis and structural stability a full causal cosmological model. The main 
result that we obtained in jit| was that only and only for s = 1/2 we obtain that the fixed point of the RG equation 
corresponds to a flat non-viscous FRW universe, i.e. a perfect fluid case. The same result is obtained under the 
standard stability analysis performed in this work, showing in this way that only for the case s = 1/2 the causal bulk 
viscous model approximates in the long time limit the dynamics of a flat perfect fluid filled FRW Universe. 


V. LIE SYMMETRIES OF THE MODEL 


In this Section we are going to study the field equations through the symmetry method. As we have indicated in 
the introduction dimensional analysis is just a manifestation of scaling symmetry, but this type of symmetry is not the 
most general form of symmetries. Therefore by studying the form of G{t) for which the equations admit symmetries, 
we hope to uncover new integrable models. But we shall see that under the assumptions made we only obtain the 
previous solutions obtained in the above sections. 

We start again with the assumption II = — xpc 2 , with x > 0. The bulk viscosity evolution equation can then be 
rewritten in the alternative form 


1 — r p 

~^~~P 


1 i_, 

~P 

a 


= 3 




(39) 


Taking the derivative with respect to the time of this equation and with the 
second order differential equation describing the time variation of the density 


use of Eq. (|l5|) we obtain the following 
of the cosmological fluid: 


p= — -Dp'-'p + AGftp 2 , 

P 

where D = > 0 and A = 12,r (7-x)(x-2 ) <- q 

q(l— r) x(l —r ) 

Equation ([f0|) is of the general form. 


P = fit, P,P), 


(40) 


(41) 


where f(t, p, p) = — Dp 1 s p + AG(t)p 2 . 

We are going now to apply all the standard procedure of Lie group analysis to this equation (see (47| for details 
and notation) 

A vector field X 


X = £{t,p)d t +rj(t, p)d p , 

is a symmetry of (^) if 


(42) 


£>ft vfp + Vtt + ( 2 T]tp £ m ) P + ivpp 2 ^) p £ppP + 


+ iv P - 2£ t - 3p£p) f-[r) t + (v P - 6) P ~ P 2 £ P \ f P = 0. (43) 

By expanding and separating ( fl3| ) with respect to powers of p we obtain the overdetermined system: 

i PP + P~^p = 0, (44) 

VP~ 2 ~ P P P~ X + V P p ~ 2 it p + 2D^ p p 1 ~ s = 0, (45) 

2?7t P - Gt + DGp 1 - 6 - 3 AGGp 2 - 2 Vt p- 1 + D (1 - s) r,p~ s = 0, (46) 

—A£Gp 2 — 2Ar]Gp + r/tt + {pp ~ 2£t) AGp 2 + Dptp 1 s = 0, (47) 

Solving (g@, we find that 

£(p, t ) = —m(l - s)t + b , r](p, t ) = mp (48) 







with the constraint 


G (2s — 1 )m 

G m( 1 — s)t — 6’ 


(49) 


and m, b are numerical constant. 

Thus we have found all the possible forms of G for which eq. (^0|) admits symmetries. There are two cases with 
respect to the values of the constant to, to = 0 which correspond to G = const, and to ^ 0 which correspond to 

G = G\ [ to .( 1 — s)t + b] 1 ~ 3 , (50) 

with G i a constant of integration. If s = | then G is a constant. For this form of G eq. ([it]) admits a single symmetry 

X = ( to ( 1 — s)t — b) dt — ( mp ) d p . (51) 

The knowledge of one symmetry X might suggest the form of a particular solution as an invariant of the operator 
X, i.e. the solution of 


dt dp 

€(t>p) v(t,p) 


(52) 


This particular solution is known as an invariant solution (generalization of similarity solution). In this case 


1 

P = p 0 t 


(53) 


where for simplicity we have taken 5 = 0, and po is a constant of integration. 

We can apply a pedestrian method to try to obtain the same results. In this way, taking into account dimensional 
considerations, from the eq. ( |40| ) we obtain the following relationships between density, time and gravitational 
constant: 


a *(1 — s)p 1 s t ^ 1, AGpt 2 1. (54) 

This last relationship is also known as the relation for inertia obtained by Sciama. From these relationship we 
obtain 

p ~ A~ 1 G ^[ 1 [to( 1 — s)t + b] 1 ~ a t ~ 2 » A~ 1 G ^ 1 [to( 1 — s)t] _TTr “ . (55) 

We can see that this result verifies the relation 

a -1 (l — s)p 1 ~ s t ^ 1. (56) 

Therefore the conditions II = — ypc 2 together with the equation © are very restrictive. Hence we have shown 
that under these assumptions there are no another possible solutions for the field equations. 


VI. CONCLUSIONS 

In the present paper we have studied a causal bulk viscous cosmological model, with bulk viscosity coefficient 
proportional to the energy density of the cosmological fluid. This hypothesis is justified by the dimensional analysis of 
the model. We have also assumed that the cosmological and gravitational constants are functions of time. With the 
help these assumptions the general solution of the gravitational equations can be obtained in an exact form, leading 
to a power-law time dependence of the physical parameters on the cosmological time. The unicity of the solution is 
also proved by the investigation of the Lie group symmetries of the basic equation describing the time variation of 
the mass density of the Universe. 

In the considered model the evolution of the Universe starts from a singular state, with the energy density, bulk 
viscosity coefficient and cosmological constant tending to infinity. At the initial moment t = 0 the relaxation time is 
r = 0. Generally the behavior of the gravitational constant shows a strong dependence on the coefficient s entering 
in the equation of state of the bulk viscosity coefficient. For s = 1/2, G is a constant during the cosmological 
evolution. From the singular initial state the Universe starts to expand, with the scale factor a an s-dependent 
function. Depending on the value of the constants in the equations of state, the cosmological evolution is non- 
inflationary {q > 0 for Hq < 1) or inflationary (q < 0 for Hq > 1). Due to the proportionality between bulk viscous 
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pressure and energy density, both inflationary and non-inflationary models are thermodynamically consistent, with 
the ratio of II and p also much smaller than 1 in the inflationary era. Generally, the cosmological constant is a 
decreasing function of time. The present model is not defined for s = 1, showing that in the important limit of small 
densities an other approach is necessary. 

Bulk viscosity is expected to play an important role in the early evolution of the Universe, when also the dynamics 
of the gravitational and cosmological constants could be different. Hence the present model, despite its simplicity, 
can lead to a better understanding of the dynamics of our Universe in its first moments of existence. 
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